Flux quantization and superfluid weight in doped antiferromagnets 
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Doped antiferromagnets, described by a t-t'-J model and a suitable 1/N expansion, exhibit a 
metallic phase-modulated antiferromagnetic ground state close to half-filling. Here we demonstrate 
that the energy of latter state is an even periodic function of the external magnetic flux threading the 
square lattice in an Aharonov-Bohm geometry. The period is equal to the flux quantum $o = 2-whc/q 
entering the Peierls phase factor of the hopping matrix elements. Thus flux quantization and a 
concomitant finite value of superfluid weight D a occur along with metallic antiferromagnetism. We 
argue that in the context of the present effective model, whereby carriers are treated as hard- 
core bosons, the charge q in the associated flux quantum might be set equal to 2e. Finally, the 
superconducting transition temperature T c is related to D s linearly, in accordance to the generic 
Kosterlitz-Thouless type of transition in a two-dimensional system, signaling the coherence of the 
phase fluctuations of the condensate. The calculated dependence of T c on hole concentration is 
qualitatively similar to that observed in the high-temperature superconducting cuprates. 

PACS numbers: 71.27.+a, 74.20.Mn 



I. INTRODUCTION 

During the last few years there has been some evi- 
dence that mobile holes in doped antiferromagnets, such 
as the high-T c superconducting copper-oxide layers [0, 
behave much like hard-core bosons. This transmutation 
of statistics, from bare fermionic holes to bosonic va- 
cancy quasiparticles, should be understood as an "emer- 
gent phenomenon" due to the reduced dimensionality 
and the presence of a strongly correlated spin back- 
ground. In the context of the simple fermionic t-J 
model, proposed by Anderson j!| to describe such sys- 
tems, the afore-mentioned evidence comes from exact- 
diagonalization studies of the ground-state energy and 
the static hole-hole correlation function on small clus- 
ters Indeed, the possibility for a hard-core boson 
behavior of the charged vacancies in doped antiferromag- 
nets, opening the way to Bose-Einstein condensation and 
the appearance of superconductivity, has been suggested 
by many authors ||[7| in the early days of the high-T c 
superconductivity research. Thoulcss in particular, 
argued that due to topological constraints, a vacancy in 
a two-dimensional torus lattice threaded by an external 
magnetic flux, must be transported twice around the ring 
in order to recover its original configuration. Hence flux 
quantization with an effective charge q — 2e may result 
from this period-doubling of the charge e bosons. 

In all the afore-mentioned works, the lack of an effec- 
tive model for doped antiferromagnets expressed in terms 
of hard-core bosons has prevented the systematic study 
of their flux quantization properties in conjunction with 
the optical and magnetic ones. Such a model, however, 
has been postulated from the outset by Psaltakis and Pa- 
panicolaou H and consists of a t-t'-J Hamiltonian and a 
suitable 1 /N expansion that provide a reasonably simple 
many-body calculational framework for the study of the 



relevant issues. When leading quantum-fluctuation ef- 
fects are taken into account in the context of this model, 
the generic experimental features of the optical conduc- 
tivity, the Drude weight, and the total optical weight in 
the cuprates are qualitatively reproduced. In particular, 
our theory |) 10 accounts aptly for the experimentally 
observed 0.5 eV peak of the midinfrared band JTT|-p^| and 
the mass enhancement factor of approximately equal to 
2 Furthermore, it predicts a finite limiting value 

for the optical conductivity <j(lu — > 0), at finite hole dop- 
ing, consistent with the residual far-infrared conductivity 
observed in the YP^CusOg+a; family of cuprates [ p^[ . 
Our results are also found to be consistent with relevant 
exact-diagonalization data |15[ . 

In view of the quoted evidence from optical experi- 
ments in favor of our effective model, we undertake in 
the present paper a systematic study of its flux quan- 
tization properties in order to provide a more complete 
assessment of the main electromagnetic responses. Our 
study includes results for the superfluid weight D s and 
the associated superconducting transition temperature 
T c . In particular, our explicit numerical estimates for 
the doping dependence of T c , including leading quantum- 
fluctuation effects, are found to reproduce qualitatively 
the observed trends in the cuprates |l^ , [r7[ |. 



II. EFFECTIVE MODEL 

Our effective model is described by a t-t'-J Hamil- 
tonian expressed in terms of Hubbard operators x ab = 
\a)(b\ as 
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where the index corresponds to a hole, the Greek indices 
fi,i>, . . . assume two distinct values, for a spin-up and a 
spin-down electron, and the summation convention is in- 
voked. Here J is the antiferromagnetic spin-exchange in- 
teraction between nearest-neighbor sites on a square 
lattice endorsed with periodic boundary conditions and 
a total number of sites A = A x x A y , where A^ = A y . 
For the hopping matrix elements Uj we assume 

!t if i,j are nearest neighbors 
— t' if i,j are next nearest neighbors (2) 
otherwise . 

The conventions in (|^) incorporate opposite signs for t 
and t' as dictated by quantum-chemistry calculations 
p8| , pj[ for Cu-0 clusters and fits of the shape of the 
Fermi surface observed by angle-resolved photoeniission 
spectroscopy p0[ . In Ref. M we generalized the local 
constraint associated with (pi) to Xi° + Xi M = ^ where 
N is an arbitrary integer, and considered the commuta- 
tion properties of the x ah 's to be those of the generators 
of the U(3) algebra. A generalized Holstein-Primakoff 
realization for the latter algebra reads 

v 00 = N — t^*?^ v' 1 " — ft 1 * t v 
Xi iv Si Si j Xi Si Si j 

(3) 

xT = (n- sTs7) 1/2 sT, xf = zt(n - ann 1 ' 2 , 

where the are Bose operators, [£f = 5ij$^ v ■ Note 
that the local constraint, giving rise to the hard-core 
character of the bosons, has been explicitly resolved in 
(||). One can then develop a perturbation theory based 
on the 1/N expansion, restoring the relevant physical 
value N = 1 at the end of the calculation. 

In the presence of an external magnetic flux <I>, thread- 
ing the two-dimensional lattice in an Aharonov-Bohm 
torus geometry, the hopping matrix elements tij are mod- 
ified by the well-known Peierls phase factor and should 
be substituted in (|I|) according to 

Uj ~» tije iA » , with Aij = -^-(Ri - Rj) • e x . (4) 

Here Ri is the position vector for site i, e x is the unit 
vector along the x-axis encircling the flux lines, and 
$0 = Zftfrc/q is the so-called flux quantum. Conven- 
tionally, the charge q of the carriers entering <3>o is, of 
course, equal to the electronic charge e. However, the ar- 
guments of Thouless |7j quoted in the Introduction imply 
that a vacancy actually "feels" twice as much external 
flux. In the context of the present effective model this 
may be accounted for by an extra factor of two in the 
expression (|J) for the which can be readily absorbed 
in a redefinition of q as q = 2e. Evidently, this reason- 
ing does not constitute a rigorous justification for the 
assignment q — 2e in the flux quantum $ - The latter 



justification can be provided only by an ab initio deriva- 
tion of an effective Hamiltonian for the hard-core boson 
vacancies, starting from a realistic electronic model for 
the cuprates. At present such a program is out of reach. 
Hence this work will be content with the study of the flux 
quantization properties of the effective model described 
by given the flux quantum constant <j> . 

In the large-TV limit "condensation" occurs, i.e., the 
Bose operators £f , \i = 1,2, become classical commuting 
fields. For uniform density states these complex number 
amplitudes may then be parametrized as 

£ = ^ e cos(^y^e-^/\ 

(5) 

= sin (|) 6 *./V*/a , 

where n e is the average electronic density, the angles Qi 
and <pi determine the local spin direction, while the re- 
maining parameter ipi determines the local phase of the 
condensate. As shown in Ref. g, close to half- filling 
(n e < 1) and for a sufficiently large t' , the ground state of 
(|l|) is described by a planar spin configuration (Oi = n/2) 
in which the local twist angles and phases are modulated 
according to 

4>i = Q • R» , ^ = Q' Ri , (6) 

where Q = (tt, tt) is the usual spin-modulating antifer- 
romagnetic wavevector and Q' = (tt, —tt) is an unusual 
phase-modulating wavevector. We should note here that 
the excitation spectrum above this ground state is gap- 
less ||, hence, as quoted in the Introduction, the lim- 
iting value of the optical conductivity a(u> — > 0) re- 
mains finite, at finite hole doping. However, as the half- 
filled-band limit is approached (n e — > 1) the quantity 
a(uj — > 0)/(l — n e ) becomes increasingly depressed [||. 
This trend is consistent with the ubiquitous "pseudogap" 
behavior observed in the optical and magnetic properties 
of underdoped cuprates pl[ | and provides further support 
to the relevance of the metallic phase-modulated antifer- 
romagnetic (AF) ground state under consideration. The 
question that now poses is how this ground state will 
respond to the presence of an external magnetic flux $? 

III. FLUX QUANTIZATION AND SUPERFLUID 
WEIGHT 

Following an argument by Yang [E2| we note that, in 
the presence of <I>, the reciprocal lattice is displaced from 
the origin by 27r$/(A x $o) along the x-axis. The quanti- 
zation of flux therefore depends on whether the ground- 
state energy of the system changes under this momentum 
boost. Given that the spin-exchange part of the Hamil- 
tonian (p])-(P|) does not couple directly to the magnetic 
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flux it is plausible that, at least in the large- TV limit, 
the condensate will respond in such a way as to leave its 
spin-modulating wavevector Q intact and simply adjust 
its phase-modulating wavevector Q' to a new value. In 
other words, we anticipate that in this classical (large- 
N) limit, the rigidity of the ground state against the 
intrusion of the external magnetic flux comes solely from 
the phase fluctuations of the condensate. These heuris- 
tic arguments lead us to consider the ansatz (Q) with the 
following modulating wavevectors 



Q = (7T,7r), Q'=(7r,-7r)- ( , Oj 



(7) 



where m is an arbitrary integer. Inserting (Q) into (^)- 
(§), the Hamiltonian takes the form = 

N 2 AEo(&), where Eo(Q) is the classical energy per lat- 
tice site for the value of physical interest N = 1. More ex- 
plicitly, taking carefully the infinite lattice limit (A — » oo) 
we have that 



A£ ($) - AE (<1> = 0) = 8iVn e (l - n e ) I — - m 

<l>o 



(8) 



Thus for each integer m we get an individual many- 
body energy level that depends quadratically on <f>. The 
ground-state energy is given by the lower envelope of 
these crossing energy-level parabolas and is characterized 
analytically by the condition 



4> 
$ 



< - , with m = 0,±l,±2,. 



(9) 



In Fig. [j] we depict by solid line the ground-state en- 
ergy calculated according to (||)-@, for typical values of 
the parameters e = t' /t, t/J, and the hole concentration 
(1 — n e ). We also depict by dashed lines the remnants of 
the individual crossing energy levels (||). Evidently, the 
ground-state energy (solid line) is an even periodic func- 
tion of the external magnetic flux $, with a macroscopic 
energy barrier between different flux minima, in accor- 
dance with the Byers and Young |2^] characterization of 
a superconductor. The period is equal to $o and there- 
fore the assignment q — 2e, discussed earlier on, leads 
to agreement with the observed flux quantization in the 
high-T c superconducting copper-oxide layers 
der to establish firmly the analytic result (t 
thus the heuristic arguments involved in 
also minimized numerically the classical energy AEq($) 
obtained by inserting directly (||) into (§)-(§), 




AE (^)=£ 1 +£ 2 

where 

E x = -n e (l - n e )^i y - 



(10) 



cos — cos -jj- cos ( Aij + 



sin — sin — cos I An 
2 2 13 



in - fj - 4>i + 4>j 



(id) 



(11) 

h) - 1] • 



The minimization of (|10|)-(^lj) was carried out by a re- 
laxation method. Excellent agreement with the ana- 
lytic result was obtained already for lattices with 
A = 20 x 20, and for all choices of the parameters e, 
t/J, and n e , within the range of stability of the phase- 
modulated AF ground state. A specific example of this 
agreement is evidenced in Fig. [l], where the open circles 
correspond to the numerical minimization data. 

Let us now turn our attention to the superfluid weight 
(or helicity modulus) D s given by the curvature of the 
infinite lattice limit of the ground-state energy AE(Q) at 

$ = o [E3IJ22U2T 



m 2 


~d 2 E{§)~ 







(12) 



$=0 



D s determines the ratio of the density of the superfluid 
charge carriers to their mass, and is related to the di- 
rectly measurable in-plane London penetration depth A l 
by D s = c 2 /(47re 2 A|). Quite generally, E(&) has an 1/N 

expansion of the form E($) = A^i^o (<f>) + NEi (<f>) H 

which leads via ( |l2|) to a corresponding expansion for the 

superfluid weight D s = N 2 Df ] + NDP H . Hence 

by exploiting the large- A^ limit result (|8|)-(||) we get im- 
mediately the expression for the leading term Di°\ 

= U'n e (l ~ ne) . (13) 

Our earlier arguments show that is a measure of the 
stiffness of the classical phase fluctuations of the conden- 
sate. Furthermore, ( p^ ) implies = Dq, where Dq 
is the leading term in the 1/N expansion of the Drudc 
weight D = N 2 D + NDi + • • • , studied in Ref. 0] 
using Kubo formalism for the current-current correla- 
tions. We have also verified, by a straightforward but 
lengthy calculation of Ei($) and the use of (12), that 

D^p — D\. Due to the analytic structure of the 1/N 
expansion, these results signify the term-by-term valid- 
ity of the identity D s = D. Strictly speaking, of course, 
we have checked explicitly that D s = D only up to and 
including terms Ds = D\, i.e., only up to and including 
leading quantum-fluctuation effects p6[ . This, however, 
is sufficient for most practical purposes and permits us 
to exploit our calculations of the Drude weight, in the 
present study of the superfluid weight. For instance, the 
weight D s = D including leading quantum-fluctuation 
effects, is found |l0| to increase linearly with small hole 
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concentration (1— n e ) away from the half-filled-band limit 
(n e = 1). This trend, present already in (p"3|), is a fun- 
damental characteristic of doped antiferromagnets. At 
higher doping values D s = D eventually saturate and 
then start to decrease. Note that the vanishing over- 
lap between the opposite sublattice spin states, along 
with the absence of quantum fluctuations in the large- N 
limit, leaves the direct hopping t' between same sublat- 
tice sites as the only relevant process of charge transport 
in this classical approximation. This argument makes 
plausible the independence of D s ^ from t and J seen in 
(|l3|). However, the leading quantum-fluctuation correc- 
tion D s — D\ involves already a non-trivial dependence 
on the latter couplings. 

It should be noted that when t' = 0, the present model 
reduces to the simple t-J model where in the physically 
relevant regime, i.e., close to half- filling, the uniform den- 
sity state under study becomes unstable against phase 
separation into an insulating (hole-poor) antiferromag- 
netic region and a conventional metallic (hole-rich) fer- 
romagnetic region (27). In the latter phase separated 
state no flux quantization and finite superfluid weight oc- 
curs. A vanishing superfluid weight has been suggested 
also by the Quantum Monte Carlo studies of the simple 
Hubbard model |p5fl , although the corresponding exact- 
diagonalization studies of the fermionic t-J model jL5| are 
not conclusive close to half-filling, due to the very small 
lattice sizes (e.g., 4x4) used. Indeed, the finite-size ef- 
fects in the numerical studies of the latter system are par- 
ticularly large because of the presence of phase separation 
in the ground state [^8| . Our observations here underline 
the importance of the next-nearest-neighbor hopping t' 
to the ability of the mobile holes in sustaining a uniform 
density state that displays flux quantization and a finite 
superfluid weight. In this respect it is useful to remind 
that the effective hopping parameter t' accounts for the 
large oxygen-oxygen overlap integrals present in the orig- 
inal CuO-2 planes 

We will complete our report with a discussion of the 
expected transition temperature to the charged super- 
fluid, i.e., superconducting, state under study. At a finite 
temperature T, the ratio of the thermal de Broglie wave- 
length of the charge carriers to their average distance 
is proportional to \JD S / (fc^T), where D s is the zero- 
temperature value determined by (|l2]). Hence a naive 
application of the criterion for the occurrence of Bose- 
Einstein condensation in an ideal boson gas, whereby 
the latter ratio should become of order unity, suggests 
a transition temperature T c of the form 



k R T r = AD, 



(14) 



type which, nevertheless, leads again to an expression 
of the form jl^ ) . Indeed, the -^-structure of the classi- 
cal Hamiltonian (10 ) — (P~l|) is a generalization of the two- 
dimensional XY model where the latter transition is well 
studied. In this context, it is important to note that a 
"universal" linear relation of the form ( |l4| ) has been es- 
tablished experimentally in the cuprates by Uemura et 
al. in their remarkable study of T c as a function of 
the zero-temperature value of A^ 2 cx D s . In the large- 
TV limit, the D s appearing in (14) is just equal to D s °' 



and the corresponding critical temperature should 
be interpreted as the ordering temperature for the classi- 
cal phase fluctuations of the condensate, in analogy with 
the analysis of Emery and Kivelson of the classical 
phase fluctuations of the conventional BCS order param- 
eter. The higher order terms in the l/N expansion of 
D s = D capture the effects of the quantum fluctuations 
and renormalize downwards these weights [ fTo| |, thereby 
reducing the corresponding value of T c . 

Following the prescription of Emery and Kivelson |29j , 
we have applied ([l4]) with A — 0.9; a numerical value 
extracted from the two-dimensional XY model [^o| . Us- 
ing the calculated D s = D of Ref. |ic|| , with the in- 
clusion of the leading quantum- fluctuation correction 
, we depict in Fig. |^ the superconducting tran- 
sition temperature T c as a function of the hole concentra- 
tion (1 — n e ). Evidently, the dependence of T c on (1 — n e ) 
reflects that of D s and reproduces qualitatively the ob- 
served trends in the cuprates 0,|l^l . With an estimated 
J/ttB ~ 1500 K in the latter materials jilll , the value of 
T c at optimum doping (1 — rig) = 0.44 (0.36), seen in 
the solid (dashed) line of Fig. |, is T c « 335 K (218 K). 
This predicted value of T c , signaling the coherence of the 
phase fluctuations of the condensate, should be regarded 
as an upper bound to an actual transition temperature 
because of the neglect of impurity disorder, higher-order 
quantum fluctuations, etc. From Fig. ^|we also note that 
with further hole doping T c starts to decrease while be- 
yond a critical doping value it vanishes, as the phase- 
modulated AF configuration, around which the present 
l/N expansion is carried out, becomes unstable. 



where A is a dimensionless constant of order unity. Of 
course, in the strictly two-dimensional model of continu- 
ous symmetry under study, a bona fide finite temperature 
phase transition can only be of the Kostcrlitz-Thouless 



IV. CONCLUSIONS 

In this paper, we have demonstrated that flux quanti- 
zation and a concomitant finite value of superfluid weight 
D s occur in the metallic phase-modulated AF ground 
state of the t-t'-J model (Q). The classical phase fluctu- 
ations of the condensate are shown to control the lead- 
ing term in the l/N expansion of D s . By appealing to 
the universality class of the two-dimensional XY model, 
the corresponding superconducting transition tempera- 
ture T c is related to D s linearly, via ( jl4| ) . The inclusion of 
leading quantum-fluctuation effects in D s provides then 
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a reasonable estimate for the order of magnitude and the 
doping dependence of T c in the cuprates. The latter de- 
pendence is of particular importance as it emerges from a 
consistent many-body 1/N expansion that preserves, at 
each order of perturbation theory, the local constraint, 
implied by the strong-correlation effects. These results 
support our effective description of the charge carriers in 
terms of hard-core bosons. 
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FIG. 1. Ground-state energy vs external magnetic flux, 
for e = 0.45, t/J = 1.0, and 1 - n e = 0.10. The zero-flux 
energy is subtracted off as to normalize the values. Solid 
line: the analytic result in the infinite lattice limit (A — * oo), 
according to Eqs. (^|)-(|^). Dashed lines: remnants of the 
crossing energy-level parabolas discussed in the text. Open 
circles: numerical minimization results for the ground-state 
energy on a finite lattice (A = 20 x 20), as determined by 
Eqs. (|lo|)-(|ll|). Evidently, the finite lattice numerical data 
(open circles) confirm the infinite lattice limit analytic result 
(solid line). 
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FIG. 2. Superconducting transition temperature vs hole 
concentration, for t/J = 1.0 and e = 0.45 (solid line) or 
e — 0.40 (dashed line), according to Eq. ( |l4| ) with the in- 
clusion of leading quantum-fluctuation effects. 
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